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DISCONTINUOUS FRACTAL FUNCTIONS AND FRACTAL 

HISTOPOLATION 

MICHAEL F. BARNSLEY AND P. VISWANATHAN 


Abstract. Fractal functions that produce smooth and non-smooth approximants constitute 
an advancement to classical nonrecursive methods of approximation. In both classical and 
fractal approximation methods emphasis is given for investigation of continuous approximants 
whereas much real data demand discontinuous models. This article intends to point out that 
many of the results on fractal functions in the traditional setting can be immediately extended 
to the discontinuous case. Another topic is the study of area matching properties of integrable 
fractal functions in order to introduce the concept of fractal histopolation. 


1. Introductory Remarks 

Approximation theory, which primarily focuses on the approximation of real-valued contin¬ 
uous functions by some simpler class of functions, covers a great deal of mathematical terri¬ 
tory. Beyond polynomials, trigonometric functions and splines, multitudes of approximation 
tools were developed which were often directed to a specific application domain. In the his¬ 
torical development of “classical” approximation theory, all efforts were directed towards study 
of smooth approximants. On the other hand, many experimental and real world signals rarely 
show sensation of smoothness in their traces and hence demand rough functions for an effective 
representation. 

Following Benoit Mandelbrot’s fractal vision of the universe, the notion of fractal function 
was introduced in reference [1] . Subsequently, it was observed that fractal functions can be used 
for smooth approximation as well [^, thereby supplementing various traditional approximation 
techniques. For more than a quarter century of its first pronouncement, theory of fractal func¬ 
tions has been extensively researched and has evolved beyond its mathematical framework; see, 
for instance, [i 1 [la [13 [Ml [H [251 [26]. Specific applications of fractal interpolation function 
include medicine, physics, and economics with, for instance, in the study of tumor perfusion 
electroencephalograms [ 23 , turbulence [ 3 , speech signals [ 21 ], signal processing [ 23 , stock 
market m etc. 

The fractal functions were originally introduced as continuous functions interpolating a pre¬ 
scribed set of data. Thus, the approximation methods, both fractal functions and their prece¬ 
dents, are based on continuity, and methods for discontinuous approximation are limited. How¬ 
ever, many real data requires to be modeled by discontinuous functions. For instance, discon¬ 
tinuous functions arise as solutions of partial differential equations describing different types 
of systems from classical physics. Since machines have only finite precision, any functions as 
represented by machine will be discontinuous and so is anything which involves discretization. 
Mandelbrot’s view of the price movement in competitive markets is one of the profound sources 
of discontinuity. In words of Mandelbrot .But prices on competitive markets need not be 

continuous, and they are conspicuously discontinuous. The only reason for assuming continuity is 
that many sciences tend, knowingly or not, to copy the procedures that prove successful in New¬ 
tonian physics....” Hence it appears that Mandelbrot envisaged not only a world of irregularity, 
but ultimate discontinuity. 

Deriving principle influence from these facts, the current article intends to investigate dis¬ 
continuous fractal functions and their elementary properties. These fractal functions may not 
interpolate the data set, thus breaking the inherent continuity and interpolatory nature of the 
fractal functions in the traditional setting. Discontinuous fractal functions reported recently in 
|19j provides one of the impetus for the current study. However, our main aim is to show that 
the tools developed in 1986 paper [T] are well-suited to deal with discontinuous fractal functions 
as well. Further, we convey that the algebraic structure of equations governing fractal functions 
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and their moment integrals are similar to that of usual continuous fractal interpolation function, 
but provides extra degrees of freedom. To this end, first we prove that the set of points of 
discontinuity of these fractal functions has Lebesgue measure zero; a fact that ensures Riemann 
integrability. As a consequence, functional equations for various integral transforms of these dis¬ 
continuous fractal functions can be easily derived, thus finding potential applications in various 
fields in science and engineering. For a restricted class of discontinuous fractal functions and its 
connection with Weierstrass type functional equations, the reader is invited to refer [3]. 

Interpolation “represents” a function by preserving function values at prescribed knot points. 
A closely related but different approximation method is Histopolation which preserves integrals 
of the function over the intervals of histopolation. Given a mesh A := {xo,xi,... ,xn} with 
strictly increasing knots and a histogram F = {/i, / 2 , ..., /a}, that is fi is the frequency for the 
interval [xi-i,Xi] with mesh spacing hi, i = 1,2,N, histopolation seeks to find a function / 
that satisfies the “area” matching condition: 



dx = fihi. 


There are several models that lead to histopolation problem. For instance, F may be obtained 
from a finite sample with the observed frequency fi in the class interval [xi-i,Xi) for i = 
1,2,... ,N and area matching function / may be taken as approximation of the unknown density 
function of the underlying random variable. The following fundamental problem occurring in one 
dimensional motion of a material point is another example where the problem of histopolation 
emerges quite naturally. Suppose we have to model velocity f{x) of the material point wherein 
g(x) is the position of the point at time x, which is known at specified points xq < xi < • • • < xn- 
Then we have 

g{x) = g{xo) + / f{t) dt, xq <t<x. 

J Xq 


Since 


f{x) dx = g{xi) - g{xi-i), i = 1,2, 


.N 


a representative for / can be obtained by solving the histopolation problem with values in the 
histogram F given by 

^ g{xi) - g{xi-i) 

Ji — 7 

hi 

In contrast to the vast literature available on interpolation, the researches on histopolation is 
limited. In case of smooth histopolant /, the problem of histopolation can be easily transformed 
in to a problem of interpolation as follows. If we construct an interpolant g € C^[xq,xn] with 
interpolation conditions g{xi) = + 5(^o)) * = 1)2,..., TV, where g(xo) is arbitrary, 

then f = g' solves the histopolation problem. This may be one of the reasons for the obscurity 
of histopolation as a separate problem. However, the situation is different in the case of fractal 
functions as they are not differentiable in general. At the same time, as in the case of interpo¬ 
lation, constructing rough histopolants is of practical relevance. Owing to these reasons, fractal 
histopolation deserves special attention and it is to this that the last section of the paper focus 
on. 

Overall, the current study may be viewed as an attempt to revitalize fractal functions and their 
applications and to initiate a study on fractal histopolation. 


2. Continuous Fractal Interpolation Function: Revisited 

To make the article fairly self-contained, we shall briefly evoke the notion of Fractal Interpo¬ 
lation Function and associated concepts in this section. To prepare the setting, first we need the 
following definition. 

Definition 2.1. Let (A, d) be a complete metric space. Let to > 1 be a positive integer and let 
Wi : X X lor i = 1, 2,..., to be continuous mappings. Then the collection {X;wi,W 2 , ■ • ■, w^} 
is called an Iterated Function System, IFS for short. 
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For a given IFS T = {X] wi, W 2 , ■ ■ ■, Wm} one can associate a set-valued map, which is termed 
collage map, as follows. Let 'H{X) denote the collection of all non-empty compact subsets of X 
endowed with the Hausdorff metric 

h(A, B) = max < maxmindla, b), maxmin(i(6, a) > VA, B G TLiX). 

^ ^ I aeA hGB ^ ' beB aeA ^ ^ ' 

Define W :H{X)-¥ 'H(X) by 

W(B) = VJZiwm, 

where Wi{B) = {wi{b) : b G B}. 

Definition 2.2. A nonempty compact subset A of A is called an attractor of an IFS J- = 
{X;wi,W 2 , ■■.,w m } if 

(1) A is a fixed point of W, that is W(A) = A 

(2) there exists an open set U C X such that A C U and limfc_>oo h(W^{B), A) = 0 for all 
B G HiU), where is the A:-fold autocomposition of W. 

The largest open set U for which (2) holds is called the basin of attraction for the attractor A 
of the IFS J-. The attractor A is also referred to as a fractal or self-referential set owing to the 
fact that A is a union of transformed copies of itself. 


If the IFS T is contractive (hyperbolic), that is each map Wi is a contraction map, then the 
existence of a unique attractor is ensured by the Banach fixed point theorem and in this case 
the basin of attraction is X. 

In what follows, the question of how to obtain continuous functions whose graphs are fractals in 
the above sense is readdressed. 

Let {[xi,yi) : i = 0,1,..., A} denote the cartesian coordinates of a finite set of points with 
increasing abscissae in the Euclidean plane Let I denote the closed bounded interval [j;o, xn] 
and li = [xi-i,Xi] for i = 1, 2,..., iV. Suppose Li : I ^ li he contraction homeomorphism such 
that 


Liixo) — T^—i, Liixj^^ — Xi^ 

\Li{x) — Li{x*)\ < li\x — T*! W x,x* G I, k G (0,1). 


Note that {/; Li, i = 1, 2,..., N} is a hyperbolic IFS with unique attractor /. Further, assume 
that Fi : J X R —>• R. is continuous map satisfying 

( Fi{xo,yo) = yi-i, Fi{xN,yN) = yi, 

I \Fi{x,y) - Fi{x,y*)\ < Si\y - y*\ y y,y* G R, Si G (0,1). 

Now define functions Wi:/xR—>-/xRby 

Wiix,y) = {Li{x),Fi{x,y)). 


( 2 . 2 ) 


The following is a fundamental theorem that leads to the definition of Fractal Interpolation 
Function. 


Theorem 2.1. (]!] j The IFS {/ x R; Wi, i = 1, 2,..., Nj has a unique attractor G{f) which is 
the graph of a continuous function /:/—>■ R satisfying f{xi) = yi for i = 0,1,..., N. 

Definition 2.3. The function / that made its debut in the foregoing theorem is termed Fractal 
Interpolation Function (FIF) corresponding to the data {(a:i, yt) : i = 0,1,..., N}. 

To obtain a functional equation for /, one may proceed as follows. 

Denote by C{I) the space of all continuous functions defined on I endowed with the Chebyshev 
norm 

llfflloo := max{|g(a;)| : x G 1} 
and consider the closed (metric) subspace 

Cvo,yN{I) ~ {g e C{I) : g(xo) = yQ-,g{xN) = 2/a}- 
Define an operator, which is a form of Read-Bajraktarivic operator, T : Cy^^y,^{I) -A Cyg^yj^(I) 
{Tg){x) = F,{L~^{x),g o L~^{x)), x G h, i G {1, 2,..., A}. 
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Theorem 2.2. (^) The operator T : Cyg^y^[I) —?> Cyg^yj^{I) is a contraction with a contractivity 
factor s := maxjsi : i = 1,2,..., N} and the fixed point of T is the FIF f corresponding to the 
data {{xi,yi) : i = 0,1,..., N}. Consequently, f satisfies the functional equation 

f{L,{x)) = F,{x,f{x)), X € I, i=l,2,...,N. 

The most widely studied FIFs in theory and applications are defined by IFS with maps 


(2.3) 


Li{x) = OiX + bi 
Fi{,x,y) = aiy + qi{x). 


where \ai\ < 1 and : / —>■ R is continuous function satisfying 


^(a^o) = Vi-i - ctiyo, qi{xN) = yi - cayN- 

The subinterval end point restraints yield 

Xi ^2—1 , 

di — , bi — . 

Xn ^0 ^0 

The parameter ai is called vertical scaling factor of the map Wi and the vector a = (oi, 02 , • ■ ■, cka) € 
(—1,1)^ is refereed to as scale vector of the IFS. 11 qi, i = 1,2,... ,N are affine maps, then the 
FIF is termed affine FIF. In this case, qi{x) = qi^x + < 7 ^ 1 , where 


Vi - Ui-I UN - yo XNyi-i - Xoyi XNyo - xoyN 

QiO — 5 Qil — 

Xn — xq Xn — xq Xn — xo xn — xo 

The following special choice of qi in (12.31) is of special interest. 

(2.4) qi{x) = ho Li{x) - aib{x). 


where the height function h is, a continuous interpolant to the data and base function 6 is a 
continuous function that passes though the extreme points (a;o,?/o) and {xM,yN)- Note that 
in case of affine FIF, h is piecewise linear function with vertices at the data points {{xi,yi) : 
i = 0,1, 2,..., TV} and & is a line joining the extremities of the interpolation interval. In con¬ 
trast to traditional nonrecursive interpolants, the FIF / corresponding to (12.31) is, in general, 
nondifferentiable. For instance, we have 


Theorem 2.3. /|I4) 1 Let {(xq, yo), (a^i, J/i), • ■ ■, (xN,yN)} be an equally spaced data set in I = 
[a:o,XAr] = [0,1]. Consider the IFS defined through the maps i2.,I\} - ll[2.4\ ), where h G C^[0,1], 
\oii\ > ^ for i = 1,2,..., N, and h'{x) does not agree with y^ — yo in a nonempty open 
subinterval of I. Then the set of points at which the corresponding FIF defined by 

f{Li{x)) = h{Li{x)) + a^{f - b){x) 

is not differentiable is dense in I. 

However, by proper choices of elements in the IFS, a fractal function / G C^(I), k G NU {0} can 
be constructed and this is the content of the following theorem. 


Theorem 2.4. (]5]} Let {(xi,yi) : i = 0,1, 2,..., be a prescribed set of interpolation data 
with increasing abscissae. Consider the IFS {I x = 1,2,...,TV}, where Wi{x,y) = 

[Li{x),Fi{x,y)'), Li{x) = OiX + bi, and Fi[x,y) = aiy + qi{x). Suppose that for some integer 
k > 0, joi] < a\, and qi € C^{I) for i = 1,2,..., N. Let 


Fi,r{x,y) 


cay + 9^^ {x) 


yi,r 


<l[''\xi) 
al — ai 


5 yN,r 


a)v o:n 


1,2,...,fc. 


If Fi-i,r{xN,yN,k) = Fi^r{xo,yo,k) for i = 2,3,... ,N and r = 1,2,... ,k, then the IFS {I x 
K; IFi : z = 1, 2,..., TV} determines a FIF f G C^(I), and is the fractal function determined 
by the IFS {I x R; {Li{x), Fi^r{x,y)),i = 1, 2,..., TV} for r = 1,2,... ,k. 
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3. Fractal Functions of More General Nature 

In this section we note that by simple modifications in the construction of continuous fractal 
interpolation function revisited in the previous section, we can break continuity and/or interpo- 
latory property of the fractal function, providing more flexibility. Although the actual fractal 
function appearing in each case discussed in the sequel may be different with different properties, 
we shall steadfastly employ the same notation / at each appearance. 


Case 1: Continuous (but not interpolatory) fractal function 


For the data {{xi,yi) : i = 0,1, 2,..., TV}, consider the IFS defined by the maps 


(3.1) 


Li{x) = ttiX + hi 

Fi{x,y) = aty + qi{x), 


where the continuous functions gi : / —>■ M are chosen such that 

{ Fi{xo,yo) = yo, 

Fn{xn, Vn) = Vn, 

Fi{x]sr,yN) = Fi+i{xo,yo) = Vi, y* G [i/i - e, + e], T = 1,2 ,..., - 1. 

Here e may be interpreted as error tolerance in measurement or noise. 


Theorem 3.1. The fractal function f corresponding to the IFS defined through i3.1\} - i3.2\) is 
continuous and satisfies \ f{xi) — yi\ < e for all i = 0,1,N. 

Proof. Consider the closed metric subspace Cyg^y^{I) := {g S C{I) : gixo) = yo,g{xN) = yN}- 
of C{I). Define T : Cy„^y^{I) Cy„^yM) by 

{Tg)[x) = F.^{L~^{x),goL~^{x)) = a.^g{L~^{x))+qi{L~^{x)), x € fi = [xi_i,Xi\, i = l,2,...,N. 

It follows at once that Tg is continuous on A = [xi-i, Xi] for each i = 1,2,..., N. Bearing in 
mind that L~^{xi) = xn and L~^^{xi) = xo, for z = 1, 2,..., iV — 1 we have 

{Tg){x~) = F,{L~^{xi),g{L~^{xi)) = FfixN,g{xN)) =FfixN,yN) = fit 

{T9){xf) = F,+i{L~^.^{xi),g{L~^.^{xi)) = F,+i{xo, g{xo)) = Fi+i{xo,yo) = Vi- 

Thus, Tg is continuous at each of the internal knots, and consequently on 7 = [a;o,T 7 v]. On 
similar lines, 

{Tg){xo) = Fi{Lf^{xo),g{Lf^{xo))) = Fi{xo,g{xo)) = Fi{xo,yo) = yo 
{Fg)(xN) = F]\[[LJf^(xN),g{LJf^(xN))) = F]\[{xN,gixN)) = FN{xN,yN) = Vn, 
demonstrating that Tg is a well-defined map on Cyg^y^{I). Furthermore, for x,y G 
\{Tg){x) - (T/i)(x)| = |a^||^(L-^(x)) - h{L-^{x))\ 

— ^||oo 

Denoting |a|oo = max{|ai| : i = 1,2,. ..,iV}, the previous inequality stipulates 

\\Tg-Th\\oo < |a|oo||y-/i||oo, 

proving contractivity of T in Chebyshev norm. Hence by the Banach fixed point theorem, T has 
a unique fixed point /. Note that /(xq) = yo, f{xN) = Vn, and 

\f{xi)-yi\ = \{Tf){xi) -y^\ = \F^{L~\xi),f{L~\xi))) - y^\ = \Fi{xN,yN) - yi\ < e, 

completing the proof. □ 


Remark 3.1. Treating yo and j/at as parameters and replacing third equation in p.2[) with the 
condition Ffix^,yN) = Fi^i(xo,yo), we obtain a continuous fractal function not attached to 
any data set. 
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Case 2: Interpolatory (but not continuous) fractal function 

For a prescribed data set {{xi,yi) : i = 0,1,..., A^}, set I = [a;o,a;Ar], li = [xi-i,Xi) for 
i = 1,2,..., N — 1 and In = sjat]. Let Li : [a;o)[a^i-ij a;i) be affine maps such that 

Li{xo) = Xi-i and Li{xJ^) = Xi for z = 1, 2,..., — 1 and Ln ■ I ^ In affine map satisfying 
LNixo) = xn -1 and Ln{xn) = xn- Observe that in contrast to the case of continuous fractal 
function, here we deal with half-open subintervals with obvious modification for the last subin¬ 
terval so that each Xi belongs to a subinterval univocally. Let qi : I ^ W he bounded function 
so that the bivariate map F) : J x R —>• IR defined by Fi{x, y) = aty + qi{x) satisfy 

Fiixo,yo) = yi-i, i = 1,2,..., N; Fjv(a;jv, 2 /a) = 2/n- 

Theorem 3.2. Consider the IFS determined by the maps Li and Fi defined in the previous 
paragraph. The corresponding fractal function f is bounded and satisfies f(xi) = z/z for i = 

0,1,...,N. 

Proof. Note that the set B{I) of all real valued bounded functions defined on I endowed with 
the supremum norm is a Banach space. Consider the closed metric subspace 

^yo.yivW = {5 G B{I) : g{xo) = z/o, g(,XN) = z/a} 

of B{I) and define T : Byg^y,.^{I) Byg^y,.^{I) by 

iTg){x) = F,{L-\x),g{Lj\x))) = aig{L-^{x)) +q^{L-^{x)), x G h, i = l,2,...,N. 

ft is plain to see that Tg is a bounded function. For z = 1, 2,..., A^, Xi-i belongs univocally to 
li and using definition of Tg 

{Tg){xi-i) = F^{L~^ {xi-i), g{L~^ {xi-i))) = F,{xo,g{xo)) = Fi{xo,yo) = Vi-i- 

Further, xn belongs to In and using A^-th piece of the definition of Tg one obtains 

iTg)ixN) = FN{Ljj^{xN),g{Ljj^{xN))) = FN{xN,g{xN)) = T/v (a;A, 2 /a) = z/a- 

Therefore, Tg is well-defined and maps into By,^^y^{I). Following the proof of previous theorem 
we assert that T is a contraction, and consequently the Banach fixed point theorem ensures 
the existence of a unique fixed point /. It follows at once that f{xi) = {Tf){xi) = yi for 
z = 0,1,...,A^, delivering the promised result. □ 

Remark 3.2. Due to the lack of “join-up” condition Fi{xN, z/a) = Fi+i{xo, yo), i = 1,2,..., N— 
1, the fractal function / in the preceding theorem has a jump discontinuity at each of the internal 
knots (and hence possibly at many other points). Same is the case, even if g^, z = 1, 2,..., A^ are 
assumed to be continuous on I. Later, to derive additional properties of bounded (discontinuous) 
fractal function, we shall assume that the maps qi, i = 1,2,..., N involved in the IFS are 
Lipschitz continuous. 

Case 3: Discontinuous fractal function 

Here we consider fractal function in B{I), which is not attached to a data set. Let {xq, xi,... ,xn} 
be a partition of / = [a;o,a;A] satisfying xq < xi < ■ ■ • < xn- As in the previous case, let 
li = [xi-i,Xi) for i = 1,2,... ,N — 1 and In = [xn-i^xn]- Let Tj : [a;o,x^r) ^ [xi-i,Xi) be 
affinities such that Li{xo) = Xi-i and Li{xN—) = Xi lor i = 1,2,..., N — 1 and Ln ■ I ^ In 
affine map satisfying Ln{xq) = xn-i and Ln{xn) = xn- For z = 1, 2,..., A^, let g^ be bounded 
function. Note that we do not require any additional conditions such as join-up conditions and 
end point conditions for the bivariate maps Fi{x,y) = aty -I- qi{x). 

Theorem 3.3. The fractal function f corresponding to the IFS defined via the maps Li and Fi 
of the previous paragraph is bounded and satisfies the self-referential equation 

f{Li{x)) = aj{x) -f q^{x). 

Proof. Define a map T : B{T) -Y B{T) by 

{Tg){x) = F^{L~^{x),g{L~^{x))) = aig{L~^{x)) +q^{L~^{x)), x G h, i = l,2,...,N. 
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Since g and m are in the linear space B[I) , it follows readily that Tg G B{I) and T is well-defined. 
As noted previously, for g,h G B{T) 

\iTg){x) - iTh){x)\ = \ai\\g{L-\x)) - h{L-^{x))\ 

< \ai\\\g - /i||oo 

^ |ci|oo II 5 ^llooi 

and hence 

\\Tg-Th\\^ < |a|oo||5 -/i||oo- 

Therefore, T is a contraction and its fixed point / enjoys the self-referential equation 

f{Li{x)) = aif{x) + qi{x), 

completing the proof. □ 


Remark 3.3. It is worth to note that instead of constant scaling factors, we may employ 
scaling functions Ofi : / —> R with suitable conditions to provide fractal functions with more 
flexibility. To to the least, we may assume that each ai is bounded and ||q!||oo := inax{||ai||oo : 
7 = 1,2,...,1V}<1 

Remark 3.4. Since fractal function / appearing in Theorem l3.3l is bounded, it belongs to C°°{I) 
and also to £^(J) for 1 < p < 00 . In particular, / is Lebesgue integrable. 


Remark 3.5. As pointed out in reference [T] for a continuous function, given / G B{I) one may 
choose b G B{I) and consider qi{x) := f[Li{x)) — aib{x). The corresponding fractal function, 
which is termed a-fractal function, denoted by /“ provides the self-referential analogue of /. In 
case b depends linearly on /, then the correspondence f ^ provides a linear operator on 
B{I), extending the notion of a-fractal operator (see [18]) to the space B{I). 


Remark 3.6. On lines similar to Theorem EH for a prescribed data set {{xi,yi) : i = 
0,1,...,A^} by proper choices of a^ and we can construct a fractal function / G B{I) 
(not necessarily continuous) that approximates the data in the sense that \f{xi) — yi\ < e 
for z = 0,1,..., A^. For instance, Lipschitz map qi can be taken as affinities qi{x) = qiox -I- qn, 
i = 1,2, ...,N and the coefficients thereof can be selected so that the function values at the 
knots 


f{xo) 


qipxo + qii 

1 — ai 


f{xN) 


qNpXN + gJVl 

1 — ajv 


f{xi) = ai+if{xo) + qi+^,pXo -I- gi+ 1 , 1 , z = 1, 2,..., A^ - 1 


are close enough to yi, i = 0,1,..., N. 


Example 3.1. We now illustrate previous results by constructing examples of continuous and 
discontinuous fractal functions which interpolate or approximate the set of data {(0,0), (0.5, 0.5), 
(1,0)}. Fig. |T] shows the affine fractal function in the standard setting (i.e., continuous and 
interpolatory) with scale vector a = (0.75, 0.75). Note that the graph of the FIF has Minkowski 
dimension D ^ 1.585 obtained as the unique solution of (see [I]) 

N 

= 1- 

i=l 

Fig. [2] represents continuous affine fractal approximants, where the scale vector is taken to be 
a = (0.5,0.5) and coefficients appearing in the affinities qi{x) = qioX + qn are chosen such that 
\f{xi) — yi\ < € with e = 0.1,0.05, and 0.005. Fig. Ejdisplays discontinuous fractal interpolation 
function corresponding to the data set wherein scaling vector is a = (0.5,0.5) and the coefficient 
Cl appearing in the affine map is chosen (at random) as i. In Fig. EJ we break both continuity 
and interpolatory conditions inherent in a traditional affine FIF. The coefficients of affinities 
are randomly chosen so that the graph passes close to the given data, that is, we construct 
discontinuous fractal function / satisfying \ f{xi) — yi\ < 0.1. 
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Figure 1. Graph of a continuous affine fractal interpolation function. 



Figure 2. Graphs of continuous affine fractal (approximation) functions. 
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Figure 3. Graph of a discontinuous affine fractal interpolation function (The 
vertical lines display discontinuities). 


4. Some Properties of Discontinuous Fractal Functions 

Recall that a continuous fractal interpolation function can be obtained either as a fixed point 
of operators defined on suitable function spaces or as attractors of IFSs. In the previous section, 
we obtained discontinuous fractal functions as fixed points of RB-operators. Here we first prove 
that the graph G(/) of this bounded fractal function / is related to the attractor of the IFS 
{/ X R; IFi, i = 1, 2,..., N}. Further, we prove that the set of points of discontinuities of / has 
Lebesgue measure zero. Our analysis follows closely (H, Theorem 3.5, p. 403) which proves a 
similar result for local fractal functions. However, we note that for the validity of the theorem 
in [7] suitable additional condition, for instance, Lipschitz continuity, is to be imposed on the 
maps Xi (which takes the role of qt in the present setting) that seems to be missing. 

As a prelude, let us recall a pair of definitions. 

If / : —>■ R is a function on a metric space {X, d) , then the oscillation of / on an open set U is 

ujf{U) = snpf{x) - inf /(x) = sup |/(a) - f{b)\ 

x^U “61/ a,beU 

and the oscillation of / at a point x* is defined as 

u)f{x*) = limtLij(Hj(x*)), 

where B^{x*) is the open ball at x* with radius e defined by B^{x*) = {x € X ■. d{x,x*) < e}. 
Note that / is continuous at x* if and only if ujf{x*) = 0. 

Theorem 4.1. Consider Wi :Jx]R—j-JxR defined by Wi{x, y) = [Lfix), Fi{x, y)) with Li and 
Fi as in Theorem \3.3[ Further assume that qt is Lipsehitz continuous with Lipsehitz constant 
Qi and Q = max{Qi : i = 1,2,..., N}. Then the IFS {I x R.; Wi, i = 1,2,..., N} is contractive 
with respect to a metric dg on defined by 

de{{x,y),{x',y')) = \x - x'\ F e\y - y'\. 
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Figure 4. Graph of a discontinuous affine fractal (approximation) function 
(The vertical lines display discontinuities). 


, „ . 1 - maxi<i<jv a, 

where U < 0 < - - - 

w 

the graph of f. 


Furthermore, the unique attractor of this IFS is the closure of 


Proof. Let {x,y), {x' ,y') G I xR. We have 

de{Wi{x,y),Wi{x',y')) = de(^{Li{x),aiy + Qiix)), {Li{x'),a^y' + gj(a;'))^ 

= \Li{x) - Li(x)\ + 0\a^y + qi{x) - {a^y' + qi{x')) \ 

< ai\x- x'\ + 6»[|aj||y - y'\ + \qi{x) - qi{x')\] 

= (oi + 0Qi)\x - x'\ + 0\a^\\y - y'\ 

< max{ai + 0Qi, |ai|}[|a: - x'\ + 9\y - y'\\. 

For 6 as mentioned in the statement of the theorem, it follows that K := m.ax{ai + 9Qi, \ai\} < 1 
and consequently that Wi, i = 1,2,..., N are contraction maps. That is, the IFS {/ x R; Wi, f = 
1,2,..., iV} is contractive, and hence it has a unique attractor, say A. Let G{f) := {(a;, /(x)) : 
X G 1} denote the graph of / and H(/ x R) denote the space of all nonempty compact subsets 
of / X R endowed with Hausdorff metric. Consider the Hutchinson operator W : 'H(/ x R) —>■ 
%[I X R) defined by 

W{B) = u(IiW,(B) 

We have to show that W(G(f)) = G{f). Let us extend W to 2^^®, the powerset of / x R. With a 
slight abuse of notation we shall denote the extension also by W. Consider {x,y) G G{f). Then 
there exists a sequence of points x„ G / such that lim„^.oo(xn,/(x„)) = {x,y). Since I = 
and li n = 0 for all i ^ j, for each fixed x„, we have Xn = Li(x(j), and hence 

(XnJiXn)) = (L*(x(j),/(L*(x(j))) = (Li(x(,), Fi(x(j,/(xj,))) 

= W.(x;,/(x:,)) G IF(G(/)) C WiMf))- 
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Consequently, {x,y) G W(G(f)) = W(G(f)), showing that 

W) c w(W))- 

Next to show that W(G(f)) C G{f), we first show that W{G) C G{f), where G := G{f) \ 
{{xi, f{xi)),i = 0,1,... ,N} . Let (a;, y) G W{G). Then {x, y) G Wi{G) for some i G {1,2,...,7V}. 
Therefore, there is an (x' ,y') such that x' G I \ {xo,xi,... ,xn}, y' = fix'), x = Liix'), and 

y = F,ix',y') = aif{x’) + q^ix') = aif{L~^ix)) +q^{L~^ix)) = f{x), 
so that (x,y) G G{f), proving IL(G') C Gif). Using this and continuity of the Hutchinson map 

Wi^) = w(g) C WiG) C 

completing the proof. □ 


The foregoing theorem shows that the graph of a bounded discontinuous fractal interpolant / 
can be approximated by the chaos game algorithm [3]. 

Remark 4.1. Note that in the case of continuous fractal function /, the graph Gif) is closed, 
and hence the attractor coincides with Gif). 

This next theorem which is just a slight variant of the Collage theorem (see [5]) hints towards 
the choices of maps Wi so that the bounded fractal function / is close to a prescribed $ G Bil). 


Theorem 4.2. Let the mappings Wi, i = 1,2,..., N in the IFS used to generate the bounded 
discontinuous fractal function f (Cf. Theorem \3.3\) he chosen such that ||$ — T$||oo < c, where 
T is the RB-operator whose fixed point yields f. Then 

ll$-/l|oo < ^ 

1 - |a|oo 

Proof. By the Banach fixed point theorem, the fractal function / = limm_>oo T'^g, where g G 
Bil) is arbitrary. Therefore we see that 

ll<i>-/l|oo = ||‘i>- lim r™$||^= lim 
m—^oo m—¥oo 

m 

= lim II 

m—¥oo 11 * ^ 

i=l 

m 

< lim 

m—>oo * ^ 

2 = 1 

m 

= lim V||T*-i($-r$)|| 

m—>oo * ^ 

2 = 1 

m 

< hm 

m ,—^oo < ^ 


2=1 


< 


1- a 


providing the assertion. 


□ 


Theorem 4.3. The set of points of discontinuities for the function f (Cf. Theorem \3.3}) is a 
Lebesgue null set. In particular, f is Riemann integrable. 


Proof. Since / is obtained by the application of Banach fixed point theorem on the Read- 
Bajraktarevic operator T, f = lim„_j.oo ffn, where = Tg„_i, and go € B)!) is arbitrary. We 
choose go = 1: the constant function on J = [a;o,a:Ar]. Note that 

gn{Liix)) = aiPn-iix) + qiix), x G I, i=l,2,...,N. 

The function gi may have finite jump discontinuities at the internal knot points xi,X 2 , ■. ■, xn-i. 
In general, g„ may have finite jump discontinuities at the internal knot points xi, X 2 ,. ■., xn-i 
of the partition and possibly at the image points Li^ o Li^ o ... Li^_„ (xr) of the internal knots. 
Therefore by a simple arithmetic, it follows that the n-th pre-fractal function g„ has at most 
iN — l)(iV"“^ -I- 1) points of discontinuity in /. Let An denote the set of points of discontinuities 
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of Qn, then \An\ < {N — 1)(7V”“^ + 1). Further, assume A = so that |A| < Kq, that 

is A is countable. We prove that / is continuous a.t x G I \ A. Choose an e > 0. Taking 
into account that qi, i = 1,2, are Lipschitz functions with Lipschitz constants Qi and 

Q = max{Qi :i = l,2,...,7V} , from the fixed point equation 

f{L,{x)) = aif{x) + qi{x) 

it follows that 

^f{Lr{I))= sup \f[L,[x)) - f{L,[y))\ 

x,y^I 

(4.1) = sup \a^{f{x) - f{y)) + q^{x) - qt{y)\ 

x,yGl 

< \a\ooUJf{I) + Q\I\, 

where |a|oo := max{|ai| : i = 1, 2,..., N} and |/| is the length of I. 

Let n = {1,2,. ..,iV} and let 11°° be the set of all infinite sequences cr = cti(72 ... of 
elements in Q. Then (11, dn), refereed to as code space, is a compact metric space, with the 
metric do defined by do(o', w) = 2“^, where k is the least index for which ak ^ Wfc. For any 
finite code a\k = aia 2 ■ ■ ■ Ok from we obtain 

(4)) = i^cri ° ^0-2 O • ■ • O Lg-fc {!)) 

< \a\ooWf{L„^ o---oL„^{I)) +Qa^-^\Il 

which on recursion yields 

(4.2) iOf{K\,{I))=\a\tiOf{I) + Q\I\ 

||a|oo - a\ 

Consider the contractive IFS {/; Ti, i = 1,2,..., N} with attractor I. The limit limfe_>oo o 
Lo-j o ... La-f. (a) is a single point independent of a S / and the coding map tt : H°° —>■ I 

TT{a) := lim o o ... (a) 

k—^oo 

is continuous and surjective. Therefore there exists a G 11°° such that 

x = 7r(cr) =n^iL^iJl). 

For any k G N, there exists a compact interval Ik such that 

xGikC 

Set J := {Ik), where = L~^ o L~^_^ o • • • o L~^. In view of (14.21) we obtain 

k 

ujfilk) = cc/(L.|JJ)) < |a|"cc/(/) + Q|J| |, , | . 

||a| - |q;|oo| 

Since / is bounded on I, |d| < xn—Xq, and 0 as fc —>■ oo, we can choose k to be large enough 
so that each summand in the previous inequality is less than |. Consequently, ujf{Ik) < e. Since 
e is arbitrary, we deduce that iOf{x) = 0, and hence / is continuous at a: € / \ A. That / is 
Riemann integrable follows now from a standard result in analysis which states that a real-valued 
bounded function / is Riemann integrable if and only if set of points of discontinuities for / has 
Lebesgue measure zero, see for instance, [22]. □ 

Theorem 4.4. Let f be the discontinuous fractal function given in Theorem Id.,91 For each 
TO S N U {0}, the moment integral 

/m = y x^f{x) da; 

can be explicitly evaluated recursively in terms of the lower moment integrals fm-i, fm- 2 , ■ ■ ■, fo, 
the scaling factors on, and the moment 

Qm = J x^Q{x) da;, 

where the function Q : / —>■ R is defined as 

Q{x) = qio L~^{x), for x G h. 
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Proof. In view of the previous theorem it follows at once that the moment integrals are well- 
defined in the Riemann sense. With a series of self-explanatory steps we have 

N 


fm = ^ f x'^fix) da; 

i=l 

N 

= XI / x^[o‘if{L ^{x)) +q^{L^ ^(a;))] da; 

i=l 

= J {aixbi)'^ f (x) dx J x'^Q{x) dx 


N m 


i=l k =0 


fk + Qr 


which may be recast as 


Also, in particular 


fm = 


SrX (T)/fc Etl + Qr. 

Q{x) dx 


/o = ^ f{x) dx = 


1 v-^A^ 

1 


Since Oi = — -we get ^ hence on account of |ai| < 1 it follows that 


A/ -1 


Xn — Xq 


□ 


We can extend the discontinuous fractal function / supported on I to whole M by defining the 
extension to zero off I, which will also be denoted by /. Now the functional equations satisfied 
by integral transforms of these discontinuous fractal functions can be easily obtained, which 
may be of interest for variety of reasons, see also [T]. The general integral transform of a “well 
behaved” function / is defined as 

/(s) = f K{x,s)f{x) dx, 

Jm. 

where K{x,s) is a suitable function, referred to as kernel of the transformation. Using the 
functional equation for fractal function / one obtains 

/(s) = J^K{x,s)f{x) dx 

N 


= X / ^(x)) +q^{L^ ^(x)] dx 

i=l 

" f 

= '^aiai / K[Li{x),s)f{x) dx-l-Q(s), 

o — 1 


where Q{s) is the integral transform of the function Q : / —>■ K defined by Q(x) = qi[L~^{x)) 
for x € li- 

This being said, it is tempting to examine transform of a fractal function with some special 
choices of kernel functions. 

Case 1: Laplace transform 
Here K{x, s) = for x > 0. Then 

N 


/(s) = ^aiai fe V(a;) dx-f ( 5 (s) 

^ — 1 I 


I 6 

N 

^ — sb- 


= '^Oiate *"^(0^5)-k ( 3 (s). 
2=1 
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Case 2: Stieltjes transform 

Taking K{x, s) = we obtain 

^ f 1 

hs) = / -7--rT/(a;) dx + (5(s) 

^ JiS-{aiX + h) 

N 

= — +Qis). 

i—1 ^ 

Case 3: Fourier transform 

For the kernel K{x^ 5 ) = ^ where j is the square root of —1, we have 


N 

fis) = X! 

N 

= aiaie^'"'^' f{ais) + Q(s). 

i=l 

In case of uniformly spaced knot sequence we obtain an explicit expression as follows. 

Denoting A(s) = ^ for equidistant knots, i.e., for = -^ previous equation yields 


dx + Q(s) 


(4.3) /(s)=A(s)/(^) + Q(s). 

Applying Equation (14.31) recursively, we obtain 


k— 1 i 


(4.4) 


/'(»)=[nA(viT)]/(i^)+i: [ n a(^)1<j{^). 


i—1 i—0 m—1 

tO 


where the empty product Om^i A( ) = 1- Note that the previous equation extemporizes 
[Equation 3.4, [T7], p. 177]. 

We have 

N 


|A(.)l<4y; 


< |a|oo, 


and therefore 


nA{j^) 


2=1 


< |oi|^ ^ 0 as A: ^ 0. 


The previous observation in conjunction with boundedness of / asserts that the first summand 
in (HU) approaches zero as fc —>■ oo. Since 




where M is such that |Q(s)| < M for all s € K, from (14.4|) we obtain 


/(4)=E[nA( 


s 

jym-l 



5. Eractal Histopolation 

Suppose that a sequence of strictly increasing knots {xo,xi,... ,XAr} and a histogram F = 
■ ■ ■, /a}, where /i € R is the frequency for the class [xi_i, Xi), i = 1,2,... ,N are given. 
For i = 1,2,..., N, let hi := Xi — Xi-i represent the step size. In fractal histopolation, we 
match average pixel intensities with our fractal function, in contrast to matching point value as 
done with interpolation. That is, we seek for an integrable fractal function / satisfying “area” 
matching condition 

(5.1) f /(x) dx = hif^. 

J Xi — x 
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Consider the IFS defined by the maps 

(5.2) Li{x) = aiX+ hi, Fi{x,y) = aty + qi{x), i = 

where \ai\ < 1 and : / —>■ K is Lipschitz continuous map. From Theorem 13.31 and Theorem 14.31 
it follows that the corresponding fractal function is Riemann integrable and satisfies 

f(x) =aif{L~^ix)) +q~^{L~^{x)), x € h, i = l,2,...,N. 

The parameters that can be varied are scaling factors ai and functions i = 1, 2,..., The 
histopolation condition prescribed in Equation (EB necessitates 


(5.3) 


hifi = j f{x) dx 

^ Ii 

= ctitti j f{x) dx + Gi J qi{x) dx 

N 

= otiGi ^ hifi +ai qi{x) dec, i = 1, 2,..., iV. 

4 — 1 


Assume Qi ^ priori fixed maps, then in the previous equation only unknown 

is the scaling factor ai, which is obtained via 


0^2 — 


hifi - ai fj qi{x) dx 
2^i—l 


i = l,2,...,N. 


However, this solution may not be feasible, since we require that |ai| < 1 for all i = 1, 2,..., A^. 

(or a less stringent condition [X]i=i ” if work in C^{I) instead of B{I)). Thus, in 

principle, the problem demands a constrained optimization. In practice, for a quicker solution, 
we can fix scaling factors ai a priori and treat g^ as unknown functions to be determined suitably 
so that the corresponding bounded integrable fractal function / satisfies Equation (EB. 


Proposition 5.1. Let a sequence of strictly increasing knots {a:o,a;i,... ,X]\[} and a histogram 
F = {/i, / 2 ,..., /jv} be given. Consider the IFS {/ x R; Wi : i = 1, 2,..., A^} defined through 
the maps given in Equation Assume that the scaling factors are selected at random so that 

joil < 1 for i = 1,2,..., N. The corresponding fractal function solves the histopolation problem 
( Cf. Equation 15.11) ) if and only if the function qt satisfies 

/ N /" / N hifi — aiGi , hifi 

(5.4) qfix)dx = — - ^^=1 , i = l,2,...,N. 

Ji a^ 

Proof. Necessary condition follows at once from Equation (j5.3p . From the functional equation 
for / one obtains (see Theorem 14.41) 


Therefore, 


f{x) dx = 


r 

E^=l ar Jj 


1-Ef= 


qfix) da; 

-1 


/ f{x) da; = aiGi / f{x) da; + / qi{x) da; 

Jii Ji JI 


Y.^=iat Sjqfix) da; 


= aiGi — ~ ~ - \-Gi I qfix) dx. 


1 1 aiai 


Substituting the stated condition on qi in the previous equation we can deduce that Jj f(x) dx = 
hifi, completing the proof. □ 
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In what follows, we shall outline some choices for qi satisfying condition in Proposition 15.11 
For instance, taking qi as affine maps qi{x) = qiox + qn, i = 1, 2,..., iV, we obtain 


N 


(5.5) hjfi = ajUi y ] hjf^ + ai(xN - xq) 

i=l 


QM 

L 2 


{XN + Xo) + qii 


*= 1 , 2 , 


,N. 


Since G (—1,1) are chosen as parameters as in the case of affine fractal interpolation function, 
the above system consists of N linear equations in 2N unknowns qio and qn, i = 1, 2,..., TV. 
Treating also as parameters we obtain 


(5.6) 


hifi O-iOli X/i=l hifi qiQ I \ • 1 O Ar 

qii = -7-- 7ri^N + xo), z = 1,2,... ,iv. 

\^N ^0) ^ 


Next suppose that we are interested to construct a continuous fractal histopolant corresponding 
to strictly increasing knots {xq, xi,, xn} and a histogram F = {/i, / 2 , ■ ■ •, /a}- Choose j/o 
and un arbitrary. Recall from Section |3] that the fractal function corresponding to Equation 
(15.211 is continuous if map i = 1,2,...,A* satisfies 

( qiixo) = 2/o(l - Oil) 
qN{xN) = */a(1 — O-n) 

ai+iHfj + gi+i(xo) = aiUN + qi{xN), z = 1, 2, 3,..., fV - 1. 


Then the next proposition follows at once. 


Proposition 5.2. Let a sequence of strictly increasing knots {xo,xi,... ,xjv} and a histogram 
F = {/i, / 2 ,..., /a} be given. Consider the IFS {/ x R; Wi : z = 1, 2,..., N} defined through 
the maps given in Equation Assume that the scaling factors are selected at random so 

that \ai\ < 1 for i = l,2,...,iV. The corresponding fractal function is continuous and solves 
the histopolation problem [Cf. Equation US. Ill ) if the function qi satisfies system of equations 
governed by Equation {5.4^ and Equation^5.7\l. 


It naturally raised the question of solvability of the system mentioned in the foregoing Propo¬ 
sition. To this end, we make some remarks. If we take qi{x) = qiox + qn and ai G (—1,1) as 
adjustable parameters, then Proposition 15.21 provides a system of 2N + 1 linear equations with 
2N + 2 unknowns yo, pN, gzo, and qn, i = 1, 2,..., TV, choosing one of the unknowns, say z/o, 
arbitrarily, the resulting square system of linear equations may be solved. Alternatively, one 
may proceed as follows. Assume values to gi(xAr), g 2 (a^A), ■ ■ ■, gA-i(a:A), 2/o, and z/at so that 
Equation dEB specifies gi(xo) and gi(xAr) for z = 1,2,... ,iV. Let gi(xo) = A, and qi{xN) = 7 ^. 
Equation 15.41 is equivalent to 


ai(xiv-xo) 

Therefore, in this case, the problem reduces to that of solving N histopolation problems with 


qi{x) dx = (xat - Xo) 
1 ca 

boundary conditions 
(5.8) 


hifi OHQi X)i=l hifi 


1 = 1,2,..., A*. 


qi{x) dx = (xjv - Xo) 


hifi - Y.f=i hif, 


ai{xN - Xo) 


qi{xo) = l3i, qi{xN)=li, z = l,2,...,Af, 


for which one can employ methods of histopolation by traditional nonrecursive functions, see, 
for instance, m- Note that the corresponding fractal histopolant / is continuous, but non- 
differentiable in general. For a special choice of qi used in the definition of a-fractal function 
(see Remark EED, construction of continuous fractal histopolant / seems to be rather easy. To 
this end, let g G C(/) and b G C{I) be such that b ^ g, b(xo) = g(xo) and b{xj^) = g{xN). 
Consider the IFS {/ x R; Wi{x, y) = (Li{x), Fi{x, y)), z = 1, 2,..., A*}, where 


Fi{x, y) = aiy Ago Li{x) - aib{x). 
Corresponding fractal function / G C{I) satisfies 

/(x) = g{x) + ai{f - b){L-\x)). 











DISCONTINUOUS FRACTAL FUNCTIONS AND FRACTAL HISTOPOLATION 


17 


Histopolation condition in Equation (15.41) reads as 

J g{x) dx — Ui 

As mentioned earlier, one can solve + 1 histopolation problems g{x) = and b{x) = 

with boundary conditions b{xo) = g{xo) and b{xN) = g{xN)- It is worthwhile to 
mention that a fractal histospline / of continuity for the knot sequence {xq, xi,..., xn} can 
be obtained by differentiating a C^+^-continuous fractal spline interpolating the data {{xi,yi) : 
i = 0,1,..., N}, where, for instance, yo = 0 and yi = yi-i + hifi for i = \,2,... ,N. Fractal 
splines interpolating a prescribed data have received much attention in the literature, see for 
example liiiniiin]- 


b{x) da; = a. 


at 

hifi ~ Ctitti 'y ^ hifi 
i=l 


Example 5.1. Consider the knot sequence {0,1/2, 1} and a histogram F = {2, 3}. We construct 
area true approximants of the histogram F by using integrable fractal functions. Consider the 
IFS defined by the maps 

,,1 ,,11 ,,1 1,,1 3 

Li[x} = -X, L 2 {x) =2^+2’ y) = 2^ + 2^+ 4’ -^ 2 ( 0 ;,y) = ^y + 2x+-. 

Here the coefficients qio and 920 appearing in the affinities are taken (at random) as 1 and 2 re¬ 
spectively, and the other coefficients are calculated using Equation ()5.6I) . Resulting discontinuous 
fractal function given by 

( y{2x) + 2x+l if a;e[0,i), 

fix) = < 

[ i/(2x - 1)-f 4a; - I if x&[\,l], 

satisfies histopolation conditions f(x) dx = 1 and fi fix) dx = |, see Figured Assume 
that the problem demands a continuous fractal histopolant. Bearing Proposition 15.21 in mind, 
taking scale factors ai = a 2 = 0.5 and assuming the value yo of the histopolant at the end point 
Xq = 0 to be 0, we solve the linear system to obtain qiq = qn = 0, ^20 = “f and ^21 = f • 
Corresponding continuous fractal histopolant satisfying the functional equation 

( |/(2x) -I- 3x if X e [0, i), 

fix) = < 

[ i/(2x - 1) - 3x-I-4 if xe[i,l], 

is depicted in Figure [6] The C^-continuous histospline in Figure [3 is obtained by differentiating 
a cubic spline fractal function g interpolating the data set {(0, 0), ( 5 ,1), (1, |)}. For details on 
cubic spline FIF, the reader may consult [giio]. Corresponding smooth histopolant satisfies the 
self-referential equation 


{ 0.4/(2x) - 1.4616x2-h0.4x-f 2.0436 if x G [0, i), 

0.4/(2x-l)-0.3384x2-hl.676x-f 0.9404 if xG[i,l]. 

Fractal histopolants may be used to model planar data with prescribed Minkowski dimension 
which controls the selection of scaling factors. Minkowski and Hausdorff dimensions of a more 
general fractal function, for instance, bounded discontinuous fractal function, continue to remain 
as an open problem. There are many other strategies for identification of free parameters in 
fractal histopolation, and quite often the particular nature of the modeling problem states the 
type of optimization to be employed. These also deserve future investigation. 
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Figure 5. Histogram and discontinuous affine fractal histopolant. 



Figure 6. Histogram and continuous affine fractal histopolant. 
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Figure 7. Histogram and C^-continuous histospline. 
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